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Chiral Optical Absorption by a Vortex in p x ±ip y -^Wave Superconductor 
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The chiral optical absorption by a single vortex in a p^iip^-wave superconductor is studied 
theoretically. The p^iipy-wave state was recently suggested as the symmetry of the order pa- 
rameter of Sr2Ru04 superconductor. Due to the violation of time reversal symmetry, there are 
two types of vortices whose winding orientation is the same or opposite to the angular mo- 
mentum of the Cooper pair. In a real material domains with p^iip^-wave states are expected. 
However, optical absorption of circular polarized light depends only on the winding of the vortex 
and has a low energy absorption peak which results in dichroism. Dichroism occurs if super- 
conductivity is realized on a single Fermi surface sheet. However, in the case of several Fermi 
surface sheets dichroism may disappear, if the both types of carriers are present, electron-like 
and hole-like. Therefore chiral optical absorption is a possible experiment to detect the orbital 
dependent superconductivity which was suggested as the superconducting state of Sr2Ru04. 

KEYWORDS: Sr2RuC>4, orbital dependent superconductivity, p^iip^-wave superconductor, vortex, chiral optical 
absorption, dichroism, domain 



-a 

G 

O 

o 



> 

o 
o 
m 
o 

0^ 
On 
-i— > 



i 

o 
o 



X 



The study of unconventional superconductivity has be- 
come one of the most attractive problems in recent con- 
densed matter research, since various examples of this 
class have been discovered among strongly correlated 
electron systems. They include some heavy fermion com- 
pounds and the high-temperature superconductors and 
the more recently discovered S^RuOzt.EP 

It has been suggested early that Sr 2 Ru04 might be 
a spin triplet p-wave superconductor, in particular, d 
to its relation to various ferromagnetic compounds B 
Meanwhile a considerable bulk of experimental evidence 
has been collected supporting this theoretical proposal. 
For example the absence of a Hebel-Slichter peak in 
NQRo' i-and the sensitivity of Tq on non-magnetic im- 
puritiesEf clearly point towards unconventional pairing. 
Moreover, the indication of broken time reversal sym- 
metry in the superconducting phase, observed in /iSR 
measurements, gives a strong argument for the p-wave 
symmetry.tf The most decisive clue for spin-triplet p- 
wave pairing comes from the Knight shift experiment 
which shows that the spin susceptibility is not affected 
by the superconducting state.EP The most likely pairing 
state with broken of the time reversal symmetry in a 
tetragonal crystal field is given by d(k)=(k x ±ik y ) z , the 
Pziipy-wave state. 

From previous investigations, it is however know that 
dichroism for the bulk states like the p x ±ip y -wave state 
are probably small B'U Hence we will concentrate here 
on the inhomogeneous regions of the superconductor the 
vortex and the surface. For the s-wave case, dichroism in 
the chiral optieaLabsorption by a single vortex was stud- 
ied previously.EHlEil' In this paper we extend this type of 
study to the p x ±ip y -wave pairing state. 

For simplicity, we assume that the superconductor is 
basically two dimensional and has a cylindrical Fermi 
surface. We then apply the method developed for the 



study of vortices in an s-wave superconductor Let 
us start with the following Bogoliubov de-Gennes equa- 
tionM 



hoUn(r)-^-{A(r)n+-f 
-Mn(r)-^-{A*(r)n- 



- □ 



.(r)] \v n (r) = E n u n (r), 
-A*(r)]}u„(r) = E n v n (r). 

(1) 



where h =-W 2 /2m-E F , n±=d/dx±id/dy, A(r) is the 
order parameter, fcp is the Fermi wave number, E n is 
the n-th excitation energy (positive value) and E-p is the 
Fermi energy (7i=l). The wave functions u n and v n de- 
scribe the quasiparticle-hole spinor state of energy E n . 
Here the vector potential is neglected and this is valid 
for strong type-II superconductors. Though Sr2Ru04 is 
not of this type, we use eq. (1) for simplicity. There is a 
relation between the solutions of E n and —E n , 

{u_ En (r), v- En (r)} «-> K n (r), u%Jr)}. (2) 

The vortex center shall be located at the origin of the 
coordinate. Then the order parameter can be written as 
A(r)=A(r)e la6 , where a=±l represents the direction of 
the vortex (positive or negative winding) and 9 is the an- 
gle of r measured relative to a;-axis. Furthermore, A(r) 
is the magnitude of the order parameter. In conventional 
superconductors the winding orientation of the vortex 
does not affect its properties apart from the direction of 
the magnetic flux. This is different in the case of the 
p^+ipy-wave superconductor, since it breaks time rever- 
sal symmetry and introduces its own winding. In the 
following we will restrict to one of the two degenerate 
states, the p^+ipy-wave states and neglect p x — ip y -wave 
paKt,which is induced near the vortex core for simplic- 
ity.y All results can be generalized through the time 
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reversal operation. 

First we address the vortex with positive winding. We 
solve eq. (1) in a finite size system of circular shape with 
radius is R. It is convenient to expand u n and v n as 

( u n (r) \ _ 1 \ - / u lnj e iw <pij(r) \ 

^=w^kr\ J ^ (3) 

where J/ is the i-th Bessel function, Zij is the j-th zero 
of J;. Substituting eq. (3) into eq. (1), we obtain, 

eijUlnj ~ i2jAj^l 2j . / U i _2j' = Ei n Ulnj, 

j' 

-tl-2jVl-2nj ~ ^2^tXjlj' u l]' = E lnVl-2nj, 
j' 

A-,^£d^ W {A W (|:4) 

+ ^ ( ^ _ 7 )A(r) ]} Wv(r) ' (4) 

where eij = (Z[j/R) 2 /2m—EF, and E\ n is the n-th exci- 
tation energy with I which is the angular momentum 
index referred to that of the spinor component u;. It 
is sufficient for our purpose to approximate the form of 
the order parameter by A(r)=Aotanh(r/£) in order to 
avoid the complications of a full self-consistence calcula- 
tion (£=z;f/Ao: coherence length; vp: Fermi velocity). 
It is important to notice that in eq. (4) ui couples with 
u;_2 only. Without the vortex, the coupling is between 
ui and vi-i due to the internal angular momentum of the 
Px+ipy-wave state. It is the phase winding of the order 
parameter around the vortex which makes the difference. 
Introducing a cut off energy luq, we can diagonalize eq. 
(4) and obtain the excitation energies as shown in Fig. 
1. The set of parameters used here corresponds to the 
quantum limit case since /cf£=4, however, the essential 
properties hold if the system is close to the classical limit 
(large fcp^). There are two kinds of bound states, one lo- 
calized around the vortex (vortex bound state) and the 
other formed near the surface at r=R (surface botuad 
state). The latter only occurs for anisotropic pairing.t3 
The vortex and surface bound states exist only for l<l. 
Due to eq. ^ the bound state with angular momentum 
index l—l contains both vortex and surface bound states 
features. 

One tool to investigate the quasiparticle states in the 
superconductor is scanning tunneling microscopy. This 
technique probes basically the local density of states, 

JV(r, E) = [\u n {r)\ 2 S(E - E n ) + \v n {r)\ 2 5(E + E n ) 

71 

(5) 

as shown in Fig. 2. Near the vortex the bounds states 
appear as distinct subgap peaks in the density of states. 
We label the bound state peaks with corresponding m 
and vi . This feature is very similar to the behavior found 
for s-wave superconductors. Different to the s-wave case, 
however, we find here that bound states yield a local 
density of states similar to the normaLmetal close to the 
surface (r=R) as shown in Fig. 2(b) £3* We observe here 
also the feature of Friedel oscillations which have a wave 



vector 2fcp. 

Next we turn to the opposite vortex (negative wind- 
ing). The Bogoliubov de-Gennes equation takes the 
form, 

etjUtnj - i^Ayy,Uy/ = Ei n Ui nj , 

j' 

-eijVinj ~ i^L&ijtj' u lj' = EnVlnj , (6) 
i' 

where u n {r) and v n (r) have been expanded as 

In this case ui couples with Naturally the surface 
bound states appear for l<0 analogous to the previous 
case. However, the vortex bound states are now re- 
stricted the angular momenta l>0, reversed to the vortex 
with negative winding. The bound state with index 1—0 
contains both a vortex and a surface bound state fea- 
tures. 

We now consider the difference in the physical prop- 
erties between the two vortices. In Fig. 4 we show 
the local density of states only near the vortex core for 
the negative winding vortex, since the winding of vortex 
does not affect it near the surface. The peak near the 
vortex core comes from zero angular momentum, since 
only Jo has a finite value at r=0. For the negative vor- 
tex wo couples with vq and they have the lowest vortex 
bound state which makes a peak at almost zero energy 
near the core. On the other hand, U2 couples with vq 
for the positive vortex case and it forms a peak near 
the core at the second lowest vortex bound state en- 
ergy (see Fig. 2). Increasing fcp£ (classical limit) the 
difference becomes small due to the reduced spacing of 
the bound state energy level. Because the vortex bound 
states are associated with specific angular momenta, we 
expect, however, to see a distinction in the chiral optical 
absorption by the vortices. The vector potential of the 
circular polarized electromagnetic radiation is described 
by A ± (r)=A q (e x ±ie y )exp[i(q-r—u>t)], where e x and e y 
are the unit vector of the x and y axes, respectively. The 
superscripts ± denotes the right and left circular polar- 
ization, respectively. We assume that q is a small wave 
vector in the z-direction. The absorption rate at 
T=0 can be expressed as 

W± = M^zZ^L^iEn+E^-u;), 

nn' 

M± n , =A q J dr[<(r)D±<,(r) - <,(r)D ± <(r)] . 

(8) 

where e is the electron charge. Since □ has the form 
of e ±l9 ( J^±ii J|), the angular momentum of u n and v n i 
must satisfy the rule l u +l v =±l (conservation of the an- 
gular momentum). The absorption of radiation causes 
the excitation of two quasiparticles. In Fig. 5 we show 
the absorption rate of the left polarized light for the vor- 
tex with positive winding. We show only three types of 
absorptions, because the weight of other processes are 
very small. The main absorption edge lies at w/Aq~1, 
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which does not depend on the vortex winding orientation 
nor the polarization of the radiation. This absorption is 
mainly caused by the excitations of a surface bound state 
and a continuum state, so that it would also appear with- 
out the vortex. Nevertheless, due to the fact that the 
surface states are a characteristic feature of anisotropic 
pairing, this low absorption edge provides evidence for 
non-s-wave pairing as the edge to the pure continuum 
excitations lies at w/2Ao— 1. 

There is a further characteristic point in the absorp- 
tion. There are two negligible contributions of low en- 
ergy. The vortex-surface state excitation is small due 
to the long distance between the vortex core and sur- 
face. In the surface-surface case, the matrix element is 
completely zero in the absence of a vortex, and, thus, de- 
pends also strongly on the distance between vortex and 
surface. Then the only low-energy contribution comes 
from the excitation of two vortex states (V-V). With the 
symmetry relation eq. (||) the selection rule can be un- 
derstood by Fig. 6. For the right circular polarized light 
there is no such contribution from the V-V absorption 
for the positive winding vortex. The absorption features 
for the vortex with negative winding is almost the same, 
if the polarization of the light is reversed as expressed 
in Fig. 6. This means that dichroism in the absorption 
is strongly connected with only the winding direction of 
the vortex relative to the polarization of the light. In the 
real material we expect that domains of the two degener- 
ate pairing states occur and that both types of domains 
can contribute to the absorption. Even if they exist we 
can observe the dichroism in the absorption by applying 
the same external field to both domains in the sample, 
since the selection of R or L polarized light depend only 
on the winding direction of the vortex—. This result is 
essentially the same as the s-wave case. Er El 

An interesting aspect occurs in connection with or- 
bital dependent supercondutivity. Various experiments 
reporting: aJarge residual density of states at low tem- 
peratures El indicate that the superconducting state of 
Sr2Ru04 is realized niedominantly on the Fermi surface 
of the 4<i xy -orbital.li§Ej- ) The other interpretation of this 
is given by a non-unitary stated If the orbital depen- 
dent superconductivity is really the case, the relevant su- 
perconducting carriers are electrons, while for the 4<i yz 
and Ad zx orbitals would also contribute hole-like carriers, 
since both electron and hole Fermi surfaces are formed 
by these two orbitals. The solutions of eq. (1) for the 
electron and hole carriers show the following symmetry 
relations: 

u^v, □ ± «-> D T , A^A*. (9) 

Equation (8) implies that the first transformation in eq. 

does not change the absorption rate. Also the second 
transformation which corresponds to p x +iPy*~^Px— Vpy 
does not affect our previous result. On the other hand, 
the third one reverses the winding of the vortex and the 
effect of polarization (R or L) such that the -condition 
for the low-energy absorption peak is reversed. El Conse- 
quently, we expect to observe the low-energy absorption 
peak for both R and L polarized light, if Ad yz and Ad zx 



orbitals contributes to the superconductivity. 

We list the result of the appearance of the low fre- 
quency peak in Table I. 



Table I. Appearance of the absorption peak at low frequency for 
the positive winding vortex. R and L represent the right and left 
circular polarization, respectively. 
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In summary, there are features in the optical absorp- 
tion which could provide evidence for unconventional 
properties of superconductivity in Sr2Ru04. One is the 
surface bound state which yields an absorption edge at 
w/Ao=l, lower than the threshold of the two contin- 
uum states excitations. The other is the dichroism in 
the absorption caused by the excitations of two vortex 
bound states. The latter can also give evidence for the 
orbital dependent superconductivity for Sr 2 Ru04. Un- 
fortunately, the energy of the vortex bound state exci- 
tation is rather low so that it would be only observable 
in the microwave range, which limits the resolution dras- 
tically. The absorption edge may, therefore, be a more 
promising experimental target to probe the unconven- 
tional superconducting state of Sr2Ru04. 
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Fig. 1. Excitation energies normalized by Ao for the negative 
vortex. S, V and C represent the surface bound state, vortex 
bound state and continuum state, respectively. Set of parameters 
are chosen as R=lOn£, fcp£=4, tt>c=2Ao. R=10tt^ corresponds 
to the level discreetness in the normal state as 5e~0.1Ao. S and 
V states close to zero angular momentum are also shown. 



Fig. 2. Local density of states for the positive vortex in an arbi- 
trary unit. A small imaginary part of 0.03Ao is added to E for 
convenience, (a) Near the vortex core, (b) Near the surface. In 
each figure r and E are normalized by £ and Ao, respectively. 
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Fig. 5. Optical absorption rate of the left circular polarized light 
by a positive vortex for the p x -Kpy-wave order parameter in an 
arbitrary unit. The frequency is normalized by Ao. A small 
imaginary part of 0.04Ao is added to u) for convenience. V-V, 
S-C and C-C represent the absorption caused by the excitations 
of vortex-vortex, surface-continuum and continuum-continuum 
states, respectively. The total absorption can be obtained by 
adding these three with the same weight. 



Fig. 3. Excitation energies normalized by Ao for the negative 
vortex. S and V states close to zero angular momentum are also 
shown. 




Fig. 4. Local density of states near the vortex core for the neg- 
ative vortex in an arbitrary unit. A small imaginary part of 
0.03Ao is added to E for convenience, r and E are normalized 
by £ and Ao, respectively. 
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Fig. 6. Transition from occupied to unoccupied bound states for 
the right (R) and left (L) circular polarized light. The order 
parameter is assumed to be p^+ipy-wave. 



